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Abstract 

This paper establishes small ball probabilities for a class of time-changed processes X o E, where 
X is a self-similar process and E is an independent continuous process, each with a certain small 
ball probability. In particular, examples of the outer process X and the time change E include an 
iterated fractional Brownian motion and the inverse of a general subordinator with infinite Levy 
measure, respectively. The small ball probabilities of such time-changed processes show power law 
decay, and the rate of decay does not depend on the small deviation order of the outer process X, 
but on the self-similarity index of X. 


1 Introduction 


Let W be a one-dimensional standard Brownian motion and let Ep be the inverse of a 
stable subordinator Dp of index /3 E (0,1), independent of W. Nane [21] established 
that the small ball probability of the time-changed Brownian motion W o Ep is given by 


( 1 ) 


P 


sup \W(Ep(t))\ < e 

0<t<l 


32T(/3) sin(/37r) 


7T 


£ 


-i) 


k -1 


- ( 2k - T 


as 


e X 0, 


where T(-) is Euler’s Gamma function and the notation f(x) ~ g(x) for two positive 
functions / and g means that lim f{x)/g{x) = 1. The result is interesting since the 
small ball probability of W o Ep shows power law decay unlike the exponential decay 
observed for the original Brownian motion W : 


( 2 ) 


— log PI sup | W (t) | < e 

m<t<i 



as e j. 0. 


Moreover, the rate of decay in (1) does not depend on the stability index f3 of the un¬ 
derlying stable subordinator Dp ; the dependence on (3 only appears as a small deviation 
constant independent of e. 

The proof of (1) provided in [21] essentially relies on the following expression for the 
Laplace transform of the random variable Ep (1) and its asymptotic behavior along the 
negative real axis (see e.g. Proposition 1(a) of [5] and Theorem 1.4 of [22]): 


( 3 ) 


E[e“ a ^ (1) ] = E p(-a) ~ 


1 

aT(l — /3) 


as a oo. 
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Here, E p(z) = Yl^Lo z n /T(n(3 + l) is the Mittag-Leffler function with parameter /?. Nane 

[21] also extended the result to a time-changed process X o Ep, where the outer process 
X is a self-similar process possessing a certain small ball probability, which particularly 
includes the case of a fractional Brownian motion. However, the exact small deviation 
constant cannot be specified unlike the situations considered in [2]; see Remark 12 for 
details of this point. 

The main motivation to analyze such time-changed processes comes from their non¬ 
standard diffusion structures. In particular, the time-changed Brownian motion W o Ep 
is non-Gaussian and non-Markovian, and is widely used to model subdiffusions, where 
particles spread more slowly than the classical Brownian particles do. Namely, the 
particles represented by the time-changed Brownian motion are trapped and immobile 
during the constant periods of the time change Ep. One interesting aspect of the time- 
changed Brownian motion is that its transition probabilities satisfy the following time- 
fractional generalization of the Fokker-Planck or forward Kolmogorov equation: 

d’tPit, x) = ^d 2 x p(t, x ), t > 0, x e K. 

Here, <9f denotes the Caputo fractional derivative operator in time of order /3 (see e.g. 

[22] ). The correspondence between the time-changed Brownian motion and the fractional 
Kolmogorov equation has been extended to those for different classes of time-changed 
processes and stochastic differential equations they drive; see e.g. [8, 9, 10, 13, 16, 17, 18]. 
The fractional Kolmogorov equations have found many applications in a wide range of 
scientific areas, including physics, [19, 27], finance [7, 15], hydrology [3], and biology [25]. 

In this paper, we establish small ball probabilities for a class of time-changed processes 
X o E, where X is a self-similar process and E is a continuous process independent of X, 
each with a certain small ball probability (Theorems 1 and 11). This largely extends the 
results in [21] in terms of both the outer process X and the time change E. Examples 
of X and E that can be handled within our framework include an iterated fractional 
Brownian motion and the inverse of a general subordinator with infinite Levy measure, 
respectively. Our strategy is to employ a version of the Tauberian theorem (Lemma 3) 
along with a general fact concerning a subordinator (Proposition 4), which is a different 
approach from what was taken in [21] to derive (1). In particular, even when E is the 
inverse of a stable subordinator, our method does not rely on the asymptotic expression 
for the Mittag-Leffler function given in (3). 

The results to be established in this paper show that the small ball probability of a 
certain time-changed process X o E has power law decay whose rate depends on the self¬ 
similarity index of the outer process X , but not on the small deviation order of X. In a 
particular case of a time-changed Brownian motion W o E with the time change E being 
the inverse of a general subordinator with infinite Levy measure, the dependence on E is 
reflected on the associated small deviation constant. We will specify that constant when 
the underlying subordinator is a Gamma subordinator or a tempered stable subordina¬ 
tor; these specific time changes have been recently investigated to analyze anomalous 
diffusions observed in various natural phenomena (see e.g. [12]). This will allow us to 
examine how the small ball probabilities for the important subclasses of time-changed 
processes vary according to the choice of the parameters defining the underlying subor- 
dinators. In particular, our result with the time change being the inverse of a tempered 
stable subordinator recovers (1) as an immediate corollary; see Remark 10 for details. 
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2 Small ball probabilities for time-changed Brownian motions 

Let E be a stochastic process in R 1 with continuous, nondecreasing paths starting at 
0. One way to construct such a process is through a subordinator. Namely, let D be 
a subordinator with Laplace exponent 0 and infinite Levy measure v, i.e. D is a one¬ 
dimensional nondecreasing Levy process with cadlag paths starting at 0 with Laplace 
transform 

poo 

(4) E[e -s£,(f) ] = where 0(s) = bs + / (1 — e~ sx ) i/(dx), s > 0, 

Jo 

with b > 0 and f^°(x A \)u{dx) < oo. The assumption that the Levy measure is infinite 
(i.e. u(0, oo) = oo) implies that 0 is an increasing function with lim s ^. oo , 0(s) = oo and 
D has strictly increasing paths with infinitely many jumps (see e.g. Theorem 21.3 of 
[24]). Let E be the inverse or first hitting time process of D] i.e. 

Eft) := infju > 0; Dfu) > £}, t > 0. 

Since D has strictly increasing paths, the process E, called an inverse subordinator , has 
continuous, nondecreasing paths starting at 0 (see e.g. Lemma 2.7 of [13]). It is known 
that E generally does not have independent or stationary increments (see Section 3 of 
[17]), which implies that even if X is a Gaussian or Levy process independent of E , the 
time-changed process X o E no longer has the same structure. Hence, existing results 
on small ball probabilities of Gaussian or Levy processes cannot be directly applied to 
find the small ball probability of X o E. 

A stochastic process X in M 1 is called a self-similar process of index H > 0 if for 
every a > 0, (X(at)) t > 0 = d (a H Xft)) t >o- Important examples of self-similar processes 
include fractional Brownian motions, iterated fractional Brownian motions, and stable 
Levy processes. Brief definitions of these processes will be provided in examples in 
Section 3. 

The following theorem largely extends Theorem 2.1 of [21] to the case when the time 
change is given by the inverse of a non-stable subordinator. 

Theorem 1. Let E be the inverse of a subordinator D with infinite Levy measure v, 
independent of a one-dimensional standard Brownian motion W. Then for all T > 0 at 
which v has no mass (i.e. u({T}) = 0), 

/ \ QO °° 1 

(5) P( sup \W(E(t))\ < e)- - u (T,oo)^2j ——-^e 2 as e f 0. 

\ 0 <t<T /IT 6 (2 k - l)' 3 

This is interpreted as P(sup 0<t<T \W(E(t))\ < e) — o(e 2 ) if v(T, oo) = 0. 

Remark 2. 1) If u(T, oo) > 0, then the small ball probability of the time-changed 
Brownian motion W o E has a power law decay. Moreover, the rate of decay of the 
small ball probability does not depend on the choice of the inverse subordinator E ; the 
dependence on E is reflected only on the constant u(T, oo). 

2) In the degenerate case when Eft) = t, clearly u = 0 and hence the small ball 
probability becomes o(e 2 ). This is because the small ball probability for the Brownian 
motion W (without a time change) has an exponential decay as in (2). 

3) If E = Eg is the inverse of a /3-stable subordinator Dp, then (5) immediately 
recovers (1). Indeed, using the explicit form of the Levy measure of Dp (see e.g. Example 
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1.3.18 of [1]), we observe that 


r°° ft t~p 

(6) 1 ' (r ’ oo) = / T t^T) x ~'~' Ax = rtwy 

When T — 1, the last expression coincides with T(/3) sin(/37r)/7r due to Euler’s reflection 
formula; consequently, the expression (5) takes the specific form given in (1). 

The proof of Theorem 1 requires some auxiliary facts to be established first. 

Lemma 3 (A version of the Tauberian theorem). Let V be a nonnegative random variable 
and let A and 6 be positive constants. Then 

E[e -ay ] ~ Aa~ 6 as a —> oo 


if and only if 

P (E < e) ~ r ^ A + ^ e e as e f 0. 

Proof. This follows from Corollary la and Theorem 4.3 of Chapter V of [26]. □ 

Proposition 4. Let E be the inverse of a subordinator D with infinite Levy measure u. 
Then for all T > 0 at which v has no mass, 

(7) F(E(T) < e) ~ u(T, oo)e as e f 0. 

This is interpreted as P (E(T) < e) = o(e) if u(T, oo) = 0. 

Proof. See Appendix. □ 

Lemma 5. Let E be the inverse of a subordinator D with Laplace exponent if. Then for 
any fixed a > 0, the Laplace transform of the function 1 H y ~K[e~ aE ^] exists and is given 
by 

(8) Ct [E[e~ aE ^}] (s) = + a )~\ s > 0. 

Proof. See Appendix. □ 

Remark 6. Lemma 5 implies that if E = Ep is the inverse of a /3-stable subordinator, 
then for a fixed a > 0, 

£,[E[e-^W]]( S )=EA. s > 0. 

Since the right hand side coincides with the Laplace transform of the function t y 
E pi-at?) (see e.g. [22]), we recover the well-known formula E[e aE > 3 T] = E j g(— at 13 ), 
which is used to derive (1) in [21]. In the proof of Theorem 1, we use (8) to guarantee 
the use of the Fubini Theorem. 


Proof of Theorem 1. By Theorem 1 of [6] (also see the proof of Theorem 2.1 of [21]), for 
all e > 0, 


(9) 


P( sup |»’(«)l<e) =-T 

M)<t<l / 7T ^ 


-1) 


k—1 


k =1 


2k — 1 


exp 


(2k - 1) V 
8^ 
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For a fixed e > 0, since E is a continuous, nondecreasing process independent of W, 
which is self-similar with index 1/2, a simple conditioning argument along with the use 
of (9) yields 


( sup \W(E(t))\ < e) = E 

pf sup bF(s) < e 

E ) 

^0 <t<T 7 

HKs<E(T) 



= E 

= un 




E 


where 


feit) := -E 
7T 


E 

L fc=i 


- 1 ) 


k -1 


2k- 1 


exp 


We also introduce the auxiliary function 


9e(t) ■= -E 


7r 


Eatrial ‘ 

L fc=1 


(2k — l) 2 n 2 E(t) 
8c 2 


(2k — l) 2 n 2 E(t) 

8^ 


Then by the Fubini Theorem for nonnegative integrands (applied to the product measure 
IP x counting measure x (e _s *dt)) and the formula (8), the Laplace transform of the 
function 1 1 —> g e (t) is given by 


A OO ^ 

c tlse](s) = - 

k =1 


E 


4 1 Ms) 

n 2k — 1 s 

k =1 


exp — 


$(s) + 


(2k - 1) 2 t r 2 E(t) 

8^ 

(2k-l) 2 n 2 ^ _1 


(«) 




8e 2 


7 r s (2k — 1) 3 7t 2 


8e 2 

< 00 , s > 0. 


This particularly implies that g e (t) < oo for (Lebesgue) almost every t > 0, but by the 
monotonicity of the function g e , we must have g e (t) < oo for all t > 0. Therefore, due 
to the Fubini Theorem, the expectation and summation in the definition of f e (t) are 
interchangeable. Thus, 


( 10 ) 



1 4 
e 2 7T 


E 

k =1 


-1) 


k -1 


2k - 1 


-E 


exp 


(2k - 1) 2 tt 2 E(T)Y 

8^ ) 


32^, (-l)^- 1 / (2k — 1) 2 7t 2 \ 

7t 3 (2k — 1) 3< ^ T y 8e 2 y ’ 


where (Pt(o) '■= aE[e~ aE ^ T ' ) ] for a > 0. By (7) along with Lemma 3, it follows that 
(Pt(o) v(T, oo) as a —> oo. Therefore, letting e j. 0 in (10) and using the dominated 
convergence theorem (which is allowed since Y2h=i V(2& — l) 3 < oo), we obtain (5). □ 

Remark 7. In the proof of (1) provided in [21], where the time change is given by the 
inverse of a stable subordinator, the asymptotic facts about the Mittag-Leffler function 
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play a significant role (see equations (2.4) and (2.5) of that paper); they are employed to 
guarantee the use of the Fubini theorem and the dominated convergence theorem. For 
the inverse of a general non-stable subordinator, however, the quantity E[e _a£:( ^] cannot 
be represented via a special function like the Mittag-Leffler function. To overcome this 
difficulty, the proof provided above employs the explicit form of the Laplace transform 
of t 14 E[e _ “ E ^] (Lemma 5) as well as a version of the Tauberian theorem (Lemma 3) 
along with a general result concerning subordinators (Proposition 4). 

Now we turn our attention to examples of time changes E which are not considered 
in [21] but can be handled by Theorem 1. This will entail small ball probabilities for 
some of the important time-changed Brownian motions representing anomalous diffu¬ 
sions observed in various fields of science. 

Let us introduce the upper incomplete Gamma function r(z, x) defined by 

poo 

r (z,x)= / e-’VMu. 

J X 


Obviously r(z, 0) coincides with the Gamma function T(z). Note that for x > 0, the 
integral defining r(z, x ) is finite even when z < 0. In particular, for (3 G (0,1) and x > 0, 
a simple application of integration by parts yields 


( 11 ) 


F(-/3,x) 


x x — T(1 — /?, x) 

P 


Example 8 (An inverse Gamma subordinator as a time change). Let E be the inverse 
of a Gamma subordinator D with parameters c, b > 0; i.e., the Laplace exponent of D 
in (4) is given by 'ijj(s) = clog(l + s/b). Then for all T > 0, using the explicit form of 
the Levy measure (see e.g. Example 1.3.22 of [1]), we obtain 

p oo poo 

u(T, oo) = / cx~ 1 e~ bx dx = c / t~ 1 e~ t dt = cT(0,bT). 

Jt JbT 

Hence, (5) with v(T, oo) replaced by cr(0,6T) yields the small ball probability of the 
time-changed Brownian motion. 

Example 9 (An inverse tempered stable subordinator as a time change). Let D be a 
tempered stable subordinator with stability index (3 G (0,1) and tempering function 
q(x), which implies that the Levy measure of D takes the form 

poo 

(12) v(dx) = x~^~ 1 q{x)dx with q(x) — / e~ Xx n(dX), 

Jo 

where // is a finite measure on (0, oo); see [23] for details. By the Fubini theorem, 

poo poo poo 

(13) z/(T, oo) = / / x~ 0 ~ 1 e~ Xx dxn(d\) = / A /3 r(- / 5, AT)/x(dA). 

Jo Jt Jo 

Note that T(—/3, A) has an alternative expression given by (11). 

Remark 10. Suppose that the tempering function q(x) in (12) is given by the simple 
exponential tilting q(x) = (3e~ Xx /T(l — j3), where A > 0 is a fixed constant. Then the 
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Laplace exponent in (4) takes the form ^(s) = (s + A)' 3 — A^, and using (11), one can 
write the constant i/(l,oo) in (13) as 


v(T, oo) 


e -\T T -p _ ^T(l -/3, AT) 

f(W) 


Letting A f 0 yields u(T, oo) = T _/3 /r(l — (3), which coincides with the constant found 
in (6) for the inverse stable subordinator; this makes sense since a tempered stable 
subordinator with the tempering factor A set to be 0 is merely a stable subordinator. 


3 Extensions 

This section establishes small ball probabilities for a large class of time-changed self¬ 
similar processes which includes the time-changed Brownian motions discussed in the 
previous section. 

Let X = (X(t))t>o be a self-similar process starting at 0 and extend X for t < 0 
using an independent copy; i.e. let X' be an independent copy of X and set X(t) := 
X'(—t) for t < 0. We call the so-defined process X = (X(t)) t ^R a two-sided process. 
Let E = (E(t)) t > o be an independent continuous process starting at 0 which is not 
necessarily nondecreasing; this implies E may take negative values. In the next theorem, 
the notation f(x) ~ g{x) means that 0 < liminf f(x)/g(x) < limsup f(x)/g(x) < oo. 
The proof employs an idea presented in the proof of Theorem 1 of [2]. 

Theorem 11. Let X be a two-sided self-similar process starting at 0 of index H > 0 
such that 

(14) — log pf sup |A"(f)| < e) & e~ T as e f 0 

MKt<l ' 

for some r > 0. Let E be a continuous process starting at 0, independent of X, such 
that 

(15) P( sup \E(t)\ < e) « e a as e f 0 

o <t<T 

for some T > 0 and a > 0. Then 

(16) P(" sup \X(E(t))\ < e) pe e a ^ H as e f 0. 

Proof. For any 6 > 0, assumption (15) is equivalent to 

P( sup \E(t)\ l/e < e) « e 6a as e | 0. 

0 <t<T 

which, by the weak order analogue of Lemma 3 (see the discussion given in Chapter V 
of [26]) with V = sup (Kt<T \E(t)\ l / e , implies that 

E [ e -a SU p 0 < t < r |i?(t)| 1 /^ _ a -6a ag a ^ OQ ' 


This is equivalent to 


(17) 


£[ e - asup °< s ,t<T l' E h)- £; (s)| 1/S j ~ a -8cr 


as a —> oo 
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due to the inequalities 

- sup | E(t) — E(s)\ < sup \E(t)\ = sup \E{t) — E{ 0)| < sup \E(t) — E(s)\. 

2 0 <s,t<T 0 <t<T 0 <t<T o <s,t<T 

Now, by assumption (14), there exist constants ci,c 2 ,e 0 > 0 such that for all e G (0,e 0 ], 

e -d6- < P f sup |x(t)| < e) < e~ C2e ~\ 

m<i< l ' 

Setting C 3 := e~ ci£ ° and C 4 := e C2£ ° , we observe that for all e > 0, 


(18) 


c 3 e~ c, ‘ ' < Pf sup | X(t)\< e) < c a e~ c " 


Let N := inf 0 <t<r E(t) and M := sup Q<t<T E(t). The assumption that -E’(O) = 0 implies 
that N < 0 and M > 0. For e > 0, using continuity of E, independence between 
(X(t)) t> 0 and (X(t)) t< o, independence between X and E, and the self-similarity of X , 
we observe that 

Pf sup \X(E(t))\ < e] = Pf sup |X(s)| < e, sup |X(s)| < e 

'0<t<T ' 'N<s< 0 0 <s<M 


= E 

= E 


Pf sup |X(s)| < e e)f( sup |X(s)| < 

/ Vi 


VAT<s<0 
pf sup |X(s)| < 


0 <s<M 


E 






E 


^o<s<i' ' '' {—N) J 

By the upper bound in (18) and the elementary inequality (x + y) p < d(p)(x p + y p ), 


x, y > 0, with p = tH, it follows that 


/h f( sup \X(E(t))\ < e) < cle- a/H ¥,\e 


-c 2 e~ T {(-N) tH +M tH } 


< cle ~ a/H E 


tH 


2 —a 


= me 


'/ h E 


-c 2 e- T (M-N) 


0 -c 2 e T sup 0< t<T \E(t)-E(s)\’’ 


(c 2 ) 


0(7 


VT,G,o{ c 2 e T )i 


where c 2 := c 2 /d(ri/), 9 := 1 /(tH) and 

(p T ,eA a ) '■= a ea E[e- asuPo ^ Tim - E{s)ll/9 }. 

Now, (17) implies that limsup^^ ipT,e,a{ a ) < 00 > an( i hence, the desired upper bound 
follows. The lower bound is obtained in a similar manner. □ 


Remark 12. 1) The rate of decay of the small ball probability of I o E in (16) does 
not depend on r appearing in (14); the information of r is reflected on the constant c 2 
introduced in the proof. 

2) Unlike Theorem 4 of [2], a simple modification of the above proof does not lead to 
a similar result concerning strong deviation orders (i.e. a result with ~ replaced by ~). 
Indeed, if we assume (instead of (14)) that 

— logPf sup |X(i)| < e) ~ ke~ T as e 0 
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for some k > 0 , then for any 8 G (0,1), we can find constants C 3 and C 4 such that 


Cse -fc(i+^ T < p( sup |x(t)| < e) < c 4 e- fc(1 - 5)e T 

MKt<l ' 


for all e > 0. This leads to 

2 2 
( 19 ) 7 -^^T, 0 ,a(cie~ T ) < e-^p^sup^ \X(E(t))\ < e) < 


where c 4 := k( 1 + 8)/d(rH), C 2 k( 1 — 8)/d(rH), and (pT,e,a(o) is as in the above 

proof. Now, if we further assume a strong deviation condition for the time change E , 
then ifT,g,a{ci) approaches a constant as a —> 00 ; however, since the constants C 3 and c 4 
depend on 8 and do not generally approach the same value as 8 —> 0 , a strong result for 
the small ball probability for X o E does not follow from (19). Note that this issue does 
not occur in the proof of Theorem 4 of [2] since the logarithmic deviation is discussed in 
that theorem. On the other hand, in Theorem 1, the explicit formula for the small ball 
probability of the Brownian motion (valid for each fixed e > 0) allowed us to establish a 
strong deviation result. 

We now consider some specific outer processes X that can be handled within the 
setting of Theorem 11. Some of the examples below show that Theorem 11 indeed 
generalizes Theorem 2.3 of [21]. 

Well-known examples of self-similar processes which have logarithmic small devia¬ 
tion orders include a fractional Brownian motion and a symmetric stable Levy process. 
Namely, if Wh denotes a fractional Brownian motion with Hurst index El G (0,1), 
i.e. Wh is a zero mean Gaussian process with covariance function E[IT#(s)Wfl-(f)] = 
(. s 2H + t 2H — |s — t\ 2H )/2, then Wh is a self-similar process of index E[ with small devi¬ 
ation order given by 


logPf sup \W H (t)\ < 

MKt<l 


c H e 


-l/H 


as e 4 , 0, 


where ch is a positive constant depending on E[. An explicit representation of the small 
deviation constant ch is found in [14]. 

On the other hand, if S a is a symmetric stable Levy process of stability index a G 
(0, 2 ], i.e. S a is a Levy process with characteristic function K[e luSa ^] = for some 

positive constant k (see e.g. [ 1 , 24]), then S a is a self-similar process of index H = 1/a 
and 

— log IP ( sup \S a (t)\ < e) ~ A a e" Q as e l 0, 

where X a > 0 is some constant; see [ 20 ] for details. 

Both of these examples satisfy condition (14) with r = l/H , and they can also be 
handled by Theorem 2.3 of [21]. However, self-similar processes with index H with 
r 7 ^ l/H also exist as the following examples show. These processes are outside the 
scope of Theorem 2.3 of [21], but Theorem 11 still applies. 

Example 13 (An iterated fractional Brownian motion as an outer process). An n- 
iterated two-sided fractional Brownian motion is the process X ^ defined by the iteration 


xW(t):=W Hl {t), X^\t):=W Hj {X^)(t)), j = 2 ,...,n, 


9 




where Wh x , • • •, Wh„ are independent two-sided fractional Brownian motions with Hurst 
indices Hi,, H n and small deviation constants ch 1 , • • •, Cff„, respectively. The process 
XW is self-similar with index := ]^[” ={ Hj. Moreover, it is established in Section 
4.2 of [2] that 

— logpf sup |A" (n )(t)| < e) ~ c n e~ Tn as e \. 0, 


where r n := 1/ ELi II"=i H j and c n is dehned iteratively by 


Ci • c-H\ , Cj . (IT 


1/Tj-i 2c Hj 


T-1 


r 7 _i J 


j = 2 ,... , 71 . 


Hence, condition (14) holds with t = r n ^ l/H^y 

Example 14 (An iterated strictly stable Levy process as an outer process). Let S ai be 
a two-sided strictly stable Levy process of index aq e (0,2], Let S Q2 be an independent 
strictly stable Levy process of index a 2 £ (0, 2] which is not a subordinator. We call the 
process X := S ai o S a2 an iterated strictly stable Levy process. It is easy to see that X 
is self-similar with index H = l/(aqa 2 ). Moreover, it is shown in Section 5 of [2] that 


— log P | sup \S ai (S a2 (t))\ <e) ~ e -a 1 a 2 /( 1 +a 2 ) £ | q 

Vo<t<i ' 


Hence, condition (14) holds with r = ai« 2 /(l T a 2 ) ^ ai« 2 = 1 / H. 

Theorem 11 also allows us to consider time changes which are given by mixtures of 
independent inverse subordinators. 

Example 15 (A mixture of independent inverse subordinators as a time change). For 
each j — 1,..., m, let Ej be the inverse of a subordinator Dj with infinite Levy measure 
Uj having no atom at T > 0 so that (7) holds for each Ej. Assume that Ej 's are 
independent and let E Y^j=\ c jEj, where Cj s are positive constants. Then it follows 
from Lemma 3 that 

771 / 771 // t n \ \ 

E[ e - a£(T) ] = Y[E[e- ac ^ E ^ T) } ~ as a 00 ’ 

j = 1 

which, again by Lemma 3, is equivalent to 

¥(E(T) < e) ~ as e | 0. 

ml \- LJ - Cj J 

\j=l o / 

Hence, Theorem 11 applies to the time change E with a = m. Moreover, with this 
specific time change, it is possible to generalize Theorems 1 to obtain the small ball 
probability of the time-changed Brownian motion with the exact small deviation constant 
specified. The proof simply combines the ideas used in the proofs of Theorems 1 and 11 
and hence is omitted. 

Note that even if each Dj is a stable subordinator, the time change E dehned in 
this example does not coincide with the inverse of a mixture of independent stable 
subordinators appearing in [9, 10, 18]. Indeed, in those papers, E is dehned to be the 
inverse of D := Xq=i where Dj 's are independent stable subordinators, which 

implies that it has the small ball probability with a = 1 due to Proposition 4. 
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Appendix 


Proof of Proposition 4- By the inverse relationship between E and D, it follows that 
P (E(T) < e) = P(-D(e) > T). Hence, we only need to verify that 


lim 

40 


P(D(e) > T) 


= u(T, 00 ) provided that z/({T}) = 0. 


Although this may be a well-known fact, for the sake of completeness of the discussion 
as well as clarification of why the assumption that u({T}) = 0 is needed, we provide a 
proof below. Note that a similar argument appears in [11], 

For a fixed real sequence {e n } with e n j. 0, let 


u n {x) 


FQD On) > a;) 

^n 


and u(x) 


00 ) 


for x > 0. Then the proof of Theorem 1.2(i) of [4] shows that the sequence of absolutely 
continuous measures u n (x)dx converges vaguely to b8o(dx) + i>(x)dx, where b > 0 is 
the drift parameter appearing in (4) and is the Dirac measure with mass at 0. This 
particularly implies that for any 0 < c < d, 


f d - 

lim / v n [x)dx 

n^oo J c 


v(x)dx. 


Now, assume that P n (T) does not converge to P(T). Then there exist a constant 
77 > 0 and a subsequence such that |P nfc (T) — P(T)| > 77 for all k. This implies that 
there exists a further subsequence {n^} such that (i) P nfcm (T) > P(T) + 77 for all m or 
(ii) P nfem (T) < P (T) — 77 for all m. If (i) holds, then since each P„ fc (•) is a decreasing 
function, for any 6 G (0, T), 


( 20 ) 


ITS 


u(x)dx = lim / v nk (x)dx > lim 8v nk (T) > 5{v(T) + 77 ). 


' T—5 


On the other hand, the assumption that z/({T}) = 0 implies that P(-) is continuous at 
x = T- hence, there exists a constant 5 E (0,T) such that v(x) — P(T) < 77/2 for all x 
with T — 5 < x < T. Thus, we have 



v(x)dx < 8{v(T) + 77 / 2 ), 


which contradicts the estimate in (20). A similar contradiction occurs if (ii) holds. 
Therefore, P n (T) must converge to P(T), which completes the proof. □ 

Proof of Lemma 5. For a fixed x > 0, since E is the inverse of D, 

F(E(t) < x) — F(D(x) >t) = 1 — P (D(x) <t), t > 0. 


Taking the Laplace transform with respect to t on both sides, we obtain 

1 1 1 _ IfrJp-s-DO)! 1 _ p -xil>(s) 

A[P (E(t) < x)](s) = - - -£t\P(D(x) E dt)](s) =- l - -^ = --- 


, s > 0 , 
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where C t [f{t)] and C t [/J,(dt)] denote the Laplace transforms of a function f(t ) and a 
measure /i(df), respectively. The right hand side of the above identity being differentiable 
with respect to x, so is the left hand side, and 

C t [P(£(t) G dx)] (s) = dx, s > 0. 

Hence, we obtain by the Fubini theorem (for nonnegative integrands) that 

C t [E[e- aE(t) ]] (s) = t^l e - xWs)+a) dx = —(?/>(s) + a)" 1 , s > 0, 

Jo s s 

which completes the proof. □ 

Acknowledgments: The author would like to thank Professor Jan Rosiiiski and Pro¬ 
fessor Xia Chen of the University of Tennessee for helpful discussions. 
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